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SECTION  I 


KINETIC  BRANCHING  OF  THE  N  *  +  02  REACTION 

\J  '  ■/ 

A  study  of  the  atmospheric  reactions  0  +  ^2  and  N  +  O2  which  yield  vibrationally  excited 
NO  molecules  has  been  initiated.  The  sources  of  NO  emission  in  the  IR  are  not  fully  understood 
and  the  detailed  kinetic  branching  of  these  reactions  is  being  studied  (References  1  and  2).  In 
particular,  the  role  of  the  N(2P]  +  02[Xvig^  reaction  in  the  production  of  vibrationally  excited 
NO  molecules  is  being  examined.  The  current  NORSE  code  assumes  that  all  of  the  exothermicity 
of  this  reaction  goes  into  vibrational  excitation  of  the  NO  product  molecule,  with  all  levels  being 
equally  populated  up  to  v  =  26.  The  purpose  of  this  report  is  to  examine  in  detail  the  kinetic 
routes  of  the  N  +  02  reaction,  including  those  which  lead  to  electronic  as  well  as  vibrationally 
excited  products.  A  preliminary  ab  initio  study  of  the  reaction  has  been  carried  out  which  suggests 
that  only  a  fraction  of  the  collisions  result  in  vibrational  excitation  in  the  product  channels. 

The  molecular  correlation  diagram  for  the  N  +  O2  reaction  is  given  in  Thble  1  for  Cs, 
and  Coov  symmetry.  This  corresponds,  respectively,  to  oblique,  perpendicular  and  linear  collisions 
of  the  N  atom  with  the  02  target  molecule.  The  adiabatic  molecular  correlations  for  C^y  and  Cs 
symmetry  are  shown  in  Figs.  1  and  2,  respectively.  These  diagrams  represent  a  more  complete 
version  of  those  given  by  Donovan  and  Husain  (Reference  3)  and  by  Schofield  (Reference  4). 

Since  the  kinetic  behavior  is  governed  by  an  average  over  all  collisional  trajectories,  we 
illustrate  the  adiabatic  correlations  in  Cs  symmetry  for  the  doublet  surfaces  in  Fig.  3,  and  for  the 
quartet  surfaces  in  Fig.  4.  In  particular,  we  illustrate  the  adiabatic  connections  from  Nf2P]  + 
02[X32Q  in  bold  lines  on  these  figures.  The  energetics  of  the  adiabatic  reaction  surfaces  arising 
from  N[2P]  +  02  and  N[2D]  +  O2  are  shown  in  Tables  2  and  3,  respectively.  For  N[2P]  +  O2, 
only  the  third  reaction  [yielding  NC^X2!!)  +  O  *  (1D)J  is  exothermic  enough  to  yield  NO  with  a 
significant  degree  of  vibrational  excitation,  provided  that  this  reaction  proceeds  along  an 
adiabatic  pathway.  However,  there  is  spectroscopic  evidence  that  curve-crossing  and 
non-adiabatic  behavior  may  occur  for  this  system.  The  ground  2A\  state  of  NO2  adiabatically 
correlates  to  the  reactants:  Nf^j  +  02[X3Z^],  whereas  the  first  excited  2Bj  state  correlates  to 
excited  state  reactants,  N[2D]  +  02[X32Q.  However,  the  2Bi  state  lies  lower  as  we  connect  to  the 
NOfX2!!)  +  0[3P]  dissociation  limit.  In  addition,  Gilmore  (Reference  5)  has  pointed  out  that  NO 
quenching  of  both  O  *  ^S)  and  O  *  [*D]  proceeds  very  rapidly,  further  indicating  curve-crossings 
in  the  doublet  surfaces. 

In  order  to  clarify  the  branching  kinetics  of  the  N  *  [2P]  +  02[X3Z^|  reaction,  we  have 
initiated  an  extensive  series  oi  calculations  of  the  potential  energy  reaction  surfaces  for  NO2  in  the 
following  symmetries:  2A',  2A",  4A',  4A".  Preliminary  Cl  calculations,  using  the  GAMESS  code 
(Reference  6),  were  carried  out  in  Cs  symmetry  assuming  a  frozen  core  of  (la’2)  (2a  )  (3a2) 
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(4a'2),  (5a'2)  (6a  2),  which  is  the  Hartree-Fock  representation  of  the  inner  Is  and  2s  electrons. 
Unfortunately,  the  characteristic  charge  density  corresponding  to  this  core  representation  differs 
significantly  for  different  spin  and  spatial  couplings  of  N  +  O2.  The  characteristic  charge  of  the 
(6a  2)  pair,  for  example,  differs  in  calculations  in  2A'  symmetry  as  compared  with  2A"  or  4A' 
symmetry.  Thus  only  a  frozen  inner  shell  (Is  electrons)  appears  to  be  possible  for  this  system. 

This  smaller  frozen  core  representation  results  in  a  great  increase  in  the  Cl  size,  from  3048  to 
28503  for  2A'  symmetry  and  from  3000  to  28125  for  2A"  symmetry.  The  calculations  must 
therefore  be  carried  out  on  a  CRAY  computer  and  the  necessary  code  modifications  for  running 
this  problem  on  the  DNA  access  machines  at  LANL  have  been  made.  The  calculations  are 
carried  out  in  Cs  symmetry  since  we  need  to  simultaneously  examine  the  N  +  O2  reaction  path  on 
nineteen  surfaces  in  2A'  and  2A"  symmetry  and  on  twelve  surfaces  in  4A'  and  4A"  symmetry. 

The  results  of  our  initial  studies  of  the  long  range  N  +  O2  behavior  are  summarized  in 
Table  4.  We  find  that  two  of  the  states  correlating  to  N[2P]  +  02[X3Z^j,  2A'X  and  2A '  IX, 
exhibit  long  range  repulsion  and  thus  should  correlate  adiabatically  to  NO[a4rT)  +  0[3P]  and 
NOpC2!!]  +  O^S]  respectively.  This  leaves  only  the  2A"  VIII  surface  correlating  to  NO[X2P]  + 
O^DJ.  For  the  quartet  surfaces  which  adiabatically  correlate  to  NO[a4IT)  +  0[3P],  no 
vibrationally  excited  ground  state  NOfX2!!]  is  predicted.  Detailed  calculations  of  these  quartet 
surfaces  are  in  progress. 

The  statistical  branching  in  the  N  +  O2  reaction  can  be  partially  analyzed  in  terms  of  these 
calculated  potential  energy  surfaces.  In  Table  2,  which  lists  the  N  *  [2P]  +  O2  reaction  surfaces, 
we  see  that  only  the  third  reaction,  which  yields  NOfX2!!]  +  0[1D],  is  exothermic  enough  to  yield 
NO  with  any  significant  degree  of  vibrational  excitation.  Thf  ^action  surfaces  for  N  *  [2D]  +  O2 
are  listed  in  Table  3.  Three  surfaces  represent  closed  (endov.iermic)  channels  for  low  energy 
collisions  and  are  non-reactive  for  air  chemistry.  Six  surfaces  produce  ground  state  NO  with 
nearly  4  eV  internal  (rovibrational)  energy  and  one  surface  produces  NO  with  less  than  2  eV 
internal  energy.  Thus,  of  the  16  potential  surfaces  governing  the  N[2D]  +  O2  and  N[2P]  +  O2 
reactions,  we  find  the  following: 

a.  3  are  non-reactive  (endothermic) 

b.  6  lead  to  electronically  excited  products  with  the  result  that  the  NO  that  is  produced  is 
not  rovibrational  excited. 

c.  7  lead  to  the  production  of  ground  state  NO  with  significant  (3-4  eV)  rovibrational 
excitation. 
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Further,  60%  of  all  collisions  involving  N[2D]  atoms  should  produce  NO  with  significant  internal 
excitation  but  only  1/6  of  all  collisions  involving  N(2P]  atoms  appear  to  lead  to  rovibrationally 
excited  NO  molecules.  Non-adiabatic  behavior  at  closer  intemuclear  separations  could  modify 
these  relative  branching  ratios  but  will  not  change  the  conclusion  that  -  1/2  of  all  collisions  of 
N[2D,  2P]  +  O2  lead  to  NO  with  little  or  no  rovibrational  excitation.  Detailed  studies  of  these 
reaction  surfaces  are  in  progress  to  further  analyze  the  branching  rates  in  this  system. 
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8'-  ' 

Figure  1.  Molecular  correlation  diagram  for  low-lying  states  of  NO2  in  symmetry. 
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Potential  energy,  eV 
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Table  1.  Molecular  correlation  diagram  for  N  +  02. 


Reactants 

Energy  (eV) 

States 

Point  Group 

N(4SU)  +  02  (X32^) 

1.401 

2,4,62  + 

Coov 

2,4,6B1 

C2v 

2,4,6A, 

Cs 

NC^)  +  O^Ag) 

2.383 

4a 

Coov 

4A2.  % 

C2v 

4A",4A' 

Cs 

N^)  +  02(b1Zg+) 

3.037 

42- 

Coov 

4a2 

C2v 

4A" 

Cs 

N(2Du)  +  ojxizj 

3.785 

2>4z+,  2>4n,  2>4a 

Coov 

i4Bi  (2),  ^B^A*  2-4A1( 

C2v 

2’4A"  (2),  2’4A'  (3) 

Cs 

N(2Du)  +  Orfa^g) 

4.766 

2s+,  2s-,  2n,  2a,  2<i>,  2r 

Coov 

2A1(2),2A2(3)>2B1(3),2B2(2), 

C2v 

2A'  (5),  2A"  (5) 

Cs 

N(2PU)  +  O^X^g) 

4.977 

2’4z-, 

Coov 

Z4Ai,  ^Az,  2,4B2 

C2v 

i4A',  Z4A"  (2) 

Cs 

N(2Du)  +  02(b1Zg+) 

5.420 

2s-,  2n,  2a 

Coov 

2A1,2A2(2)>2B1,2B2 

C2v 

2A'  (2),  2A"  (3) 

Cs 

H(X)  +  02  (C1^) 

5.499 

42  + 

Coov 

4Ai 

C2v 

4A' 

Cs 

N(4Su)  +  C>2(A'  3au) 

5.701 

2.4,6a 

Coov 

2,4,6 A 1,  2,4,6B2 

C2v 

2,4,6A',  2.A6A* 

Cs 

H(%)  +  02(A3Zu  +  ) 

5.789 

2.4,62- 

Coov 

2,4’6B2 

C2v 

2,4,6  A  " 

Cs 
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Table  1.  Molecular  correlation  diagram  for  N  +  02.  (Concluded) 


Reactants 

Energy,  (eV) 

States 

Point  Group 

N(2PU)  +  02(a1A&) 

5.958 

^A.2* 

Coov 

2Ai(2),  2 A 2,  2Bh  2B2  (2) 

C2v 

2A'  (3),  2A"  (3) 

Cs 

N(2PU)  +  02  (b12g+) 

6.612 

2z+,2n 

Coov 

2A,,  2Blt  2B2 

C2v 

2A'  (2),  2A" 

Cs 

Products 

Energy  (eV) 

States 

Point  Group 

NCK2]!)  +  0(3Pg) 

0 

2’42+,2>42-  *4A 

Coov 

i4A'(3X  2,4A"(3) 

Q 

NCX2!!)  +  O  (JDg) 

1.967 

2Z  +  ,  22-,  2T!(2),  2A,  2<D 

Coov 

2A'  (5),  2A"  (5) 

Cs 

NCK2!!)  +  O  (iSg) 

4.189 

2n 

Coov 

2A',  2A" 

Cs 

NO  (a4n)  +  O  (3Pg) 

4.765 

Z4.62  +  ;  2,4.62-  2,4.6n,  2.4, 6A 

Coov 

Z4*6A’  (3),  2’4-6A"  (3) 

Cs 

NO(A2Z  +  )  +  0(3Pg) 

5.450 

2’4z-,  2>4n 

Coov 

2’4A',  2’4A"  (2) 

Cs 

NO  (B2!!)  +  O  (3Pg) 

5.693 

2’4z\2’4x-,2’4n,2’4A 

Coov 

2’4A'  (3),  2’4A"  (3) 

Cs 

NO  (b4Z~)  +  O  (3Pg) 

6.0344 

2’4,6z + ,  2’4>6n 

Coov 

2>4  6Ai  (2),  2'4>6A" 

Cs 
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Table  2.  Adiabatic  reaction  surfaces  arising  from  N[2P]  +  O2. 


N  *  [2P]  +  02[X3S2 


N02*[2A'  X] 
N02  *  [4A'  VI] 
N02*  [2A"  VIII] 
N02‘[4A"  V] 
N02*[2A"  IX] 
N02*  [4A"  VI] 


NO  *  [a4IT]  +  O  [3P]  +  0.212  eV 
NO  *  [a4rT]  +  O  [3P]  +  0.212  eV 
NO  [X^  +  0*[1D]+  3.010  eV 
NO  *  ]a4n]  +  O  [3P]  +  0.212  eV 
NO  [X2!!]  +  O’pS]  +  0.788 eV 
NO  *  [a4n]  +  O  [3P]  +  0.212  eV 


Table  3.  Adiabatic  reaction  surfaces  arising  from  N[2D]  +  02. 


N  *  [2D]  +  02[X32  j 


N02*[2A'  II] 
N02  *  (2A'  III] 
N02*(2A'  IV] 
N02*(2A"  I] 
N02*[2A  "  II] 
N02  *  [4A'  III] 

no2  *  l4A'  rv] 
N02,[4A'  v] 
N02*[4A"  hi] 
N02*  [4A"  IV] 


NO  [ X zrn  +  O  [3P]  +  3.785  eV 

NO  pC2^  +  O  [3P]  +  3.785  eV 

NO  [X2!!]  +  O  *  [2D]  +  1.878  eV 

NO  [X2!!]  +  O  [3P]  +  3.785  eV 

NO  {X2^  +  O  [3P]  +  3.785  eV 

NO  [X2^  +  O  [3P]  +  3.785  eV 

NO  *  (a4ri]  +  O  ]3P]  -  0.980  eV 

NO  *  [a4n]  +  O  [3P]  -  0.980  eV 

NO  [X2^  +  O  [3P]  +  3.785  eV 

NO  *  [a4IT]  +  O  [3P]  -  0.980  eV 
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Table  4.  Minimum  energy  long  range  interaction  potentials  for 
N  +  O2  in  Cs  symmetry. 


Long  Range 

2A'  R  =  4.0  bohrs*  R  =  3.0  bohrs  Character 


4Su  +  X32g 

-201.4657103 

-201.4657669 

a 

-201.3423547 

-201.3428036 

a 

2Du  +  X32g 

-201.3423540 

-201.342542 

r 

-201.3423535 

-201.3415228 

r 

-201.3045858 

-201.3051402 

a 

-201.3045827 

-201.3048667 

a 

2Du  +  a3Ag 

-201.3045808 

-201.3043220 

r 

-201.3045790 

-201.3041826 

r 

-201.3045775 

-201.3040453 

r 

2PU  +  X3!' 

-201.2994959 

-201.2988189 

r 

2A 


2Du  +  X32g 


2Du  +  a!Ag 
2PU  +  X3Zg 


-201.3423546 

-201.3427958 

a 

-201.3423536 

-201.3420978 

r 

-201.3045845 

-201.3051136 

a 

-201.3045830 

-201.3048777 

a 

-201.3045810 

-201.3043588 

r 

-201.3045790 

-201.3041614 

r 

-201.3045775 

-201.3040561 

a 

-201.2994986 

-201.2999683 

a 

-201.2994921 

-201.2987603 

r 

aR  is  measured  from  N  to  the  center  of  the  O2  molecule,  energy  in  hartrees. 
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